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Abstract 



In this paper, we presented a model of pulse-coupled oscillators distributed on a bidirectional ring 
' with propagation delay. In numerical simulations based on this model, we observed phenomena of 

U: asynchrony in a certain range of delay factor a. To find the cause of these phenomena, we used a new 

probability based approach of analyzing. In this approach, the mathematical expectation of influence 
on one oscillator's phase change caused by its neighbor, which is regarded as a random factor, is cal- 
culated. By adding this expectation of influence into the firing map h(cf>) introduced by Mirolla and 
Strogatz, a probability firing map H is invented. By observing the behavior of H 's iteration from H 's 
graph, we successfully constructed a connection between the asynchrony phenomena and H 's graph. 
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O 1 I. Introduction & Review 



The phenomena of mutual synchronization widely exist in nature. For example, the synchronization of 
flashes of fireflies; the unison in which crickets and cicadae chirp; the pacemaker of the heart; and the 
synchronization of women's menstrual period who live together. 

A lot of researches have been done on this topic since last century. Mirolla and Strogatz[l] has 
studied a model of oscillators' mutual synchronization with all-to-all connection and no delay pulse 
propagation. In their model, they introduced a phase-state function / : [0, 1] — > [0, 1], which is required 
to be smooth, monotonic increasing, and concave down. With this /, the phase change of an oscillator 
receiving firing was described as 4> t+ = f^ 1 {f{4> t +e)), where e is the coupling strength. In their article, 
the initial conditions in wh ich the oscillators function will not achieve synchrony, was proved to be of 
Lebesgue measure zero, given the oscillators' phase-state function satisfies the condition mentioned 
before. 

Later, the dynamics synchronization under a threshold restriction <fr c was analyzed by Chia-Chu Chen 
[2]. With threshold included, the phase of an oscillator will not be "pulled up" by others' firing when 
its phase ^ <f> c . Chen concluded that with <fi c < |, synchronization can occur for almost all initial 
conditions. 
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In Mathar and Mattfeldt's paper[3], the same conclusion was drawn on a collection of oscillators in 
all-to-all connections, but the differentiable constraint on the phase-state function introduced in[l] was 
removed. Instead / was only required to be strictly concave and increasing. They proved the conclusion 
on two oscillators' case by analyzing the iteration of the firing map h(<f>) — f~ 1 (f(l — <fi) + e) and its fixed 
points. The conclusion was expanded the cases of a clique of N oscillators by utilizing Fubini's law. In 
the second part of the article, a propagation delay a was added into the model. A similar conclusion 
that oscillators with all the initial states except a Lebesgue measure zero set will achieve synchrony was 
reached on the model with delay. 

Similar to this paper, a case of N oscillators located on a unidirectional ring was discussed in [4]. 
In their model an oscillators pulled up to firing phase will fire, while in Mirolla & Strogatz's model[l] 
it won't. As a result of this mechanism, an "avalanche" effect exists (a chain of oscillators firing one 
immediately after another), and so helps the synchronization of the ring. 

Diaz-Guilera, Perez and Arenas[5] described the FD (firing + driving) process in the form of a vector's 
multiplying a transformation matrix M. Thus, the evolution of states of oscillators was evaluated by 
multiplying a series of matrices on the phase vector (fa, fa, . . . <Pn) T ■ By comparing M's eigenvalues 
with 1, the expected evolution result of the oscillators was classified into converging (repeller fixed 
points) and phase-locking (attractor fixed points). Also the properties of fixed points was related with 
sign of e. Similar methods were used in [6], which analyzed pattern formation on a ring of oscillators 
with couple strength e < (denoted as p in this paper)selection. 

In this paper, we first presented the model of oscillators distributed on a bidirectional ring with 
propagation delay. Based on this model, numerical simulation were made. From the results of the 
simulation, we spotted phenomena of asynchrony happening in a certain rang of propagation delaying 
factor a (see Fig 7 on page 6). In purpose of explaining the asynchrony phenomena, we divided the 
analysis into three parts: 

Step 1 For a pair of neighboring oscillators A and B, once they achieve a-synchrony, their synchroniza- 
tion is easy to restore regardless of the influence from other neighbors. (Theorem 1 on page 8) 

Step 2 For a neighborhood consisted of oscillators B,C,D on the ring, we treat D as a random factor 
which uniformly distributed in [0, 1), thus we replace D with its mathematical expectation of in- 
fluence on C. 

Step 3 In a four-element neighborhood (A, B, C, D), if A and B have achieved a-synchrony once, then 
we ignore A's effect on B using the conclusion established in step 1. Also we substitute D with its 
mathematical expectation of influence on C calculated in step 2. After dealing with A and D, we 
are able to trace the synchronization process between B and C (See Formula 8 on page 13). 

By studying the iteration behavior of the firing map, which describes phase change of B and C, we suc- 
cessfully made a connection between the asynchrony phenomena (Fig 7 on page 6) and the graph of the 
firing map (Fig 10 on page 17). 



II. The Model 



In this paper we will discuss a system of N oscillators distributed on a ring. This system includes the 
following properties: 

Phase For each oscillator i, there is a real number e [0, 1], called a phase variable, belongs to it to depict 
its phase. 

Driving As time goes by the phase fa of each oscillator increases spontaneously at a constant speed 
-^j- = S = ^> nere T i s tne cycle period. For convenience we let S = 1, by which we will also 
haveT = 1. See [1]. 

Firing When an oscillator i's phase reaches 1, it will fire a pulse to other oscillators immediately. ([1]) 
Ring Topology The firing of an oscillator will be received only by those oscillators that have connection 

to the firing one. In the bidirectional ring topology of this model, i's pulse is received by j only if 
j = i - 1 or i + 1 (modTV). See Fig 1 on the next page. 
Delay Similar to [3], a factor of reaction and propagation delay is included in this model: an oscillator fired 

upon will not react to the pulse until a time delay a passes, (i.e. i fires at to, j only reacts at to + a). 
Period of no responding The oscillator will not react if its phase is in the interval [0, 2a) i.e. j will not 

react if <j>j € [0, 2a) at t + a. 
Coupling When an oscillator j reacts to a pulse (with precondition, 0* ^ [0, 2a) at t as the pulse is 

fired at t — a), its phase is "pulled forward" for p ■ </> j + e. That is 



Here p denotes the coupling strength, and e the coupling constant. In this model we always have 

p > 0, e > 0. Note that when an oscillator's <j> is "pulled" up by others to reach 1, it will still fire. 
Phase difference The phase difference between two oscillators i and j are defined as 



With this definition, the difference between cpi = 0.9 and <j>j = 0.1 is 0.2, not 0.8. 
a-synchronization When oscillators i and j satisfy E)(0j, <f>j) ^ a, they are asynchronized, or in other 
words, achieved a-synchronization. See [3]. In addition, in a system of oscillators distributed 
on a ring, this system achieves a-synchronization when each pair of "neighbor" achieves a- 
synchronization: Mi, if j = i — 1 or i + 1 (mod N) => 4>j) ^ a. A system that has achieved 

a-synchronization is showed in Fig 3 on the following page. 
Respond only once An oscillator responds to a pulse from the same "neighbor" only once in one firing 

cycle, e.g. If in some way i "catches" up with j and fires the second time before j's firing, j does 

not respond to i's second pulse. So far no previous models has added this restriction. 



(1 + p) ■ <t>f + e 




if (1 + p) ■ <j>) + e < 1 
if (1 + p) ■ <j>y +e>l 




III. Numerical simulation 



Before presenting a theoretical analysis on the dynamics of the system, first we provide some results on 
the numerical simulation. 




Figure 1: A ring of oscillators represented by cylinders. Figure 2: A system in its synchronization process. 
The height of cylinder depicts the phase of an oscillator. 
The connection lines on the bottom indicate the neighbor- 
ing relationship. 




Figure 3: A system that achieved a-synchronization 

a. Process of synchronization 

The first experiment tracks the process of synchronization. In purpose of inspecting the course of a 
system's synchronizing, a demonstration program is written in MATLAB, which shows the dynamics 
of oscillators represented by cylinders. Fig 1, 2, ?? and 3 display four snapshots in a synchronization 
procedure. 

Also the count of a-synchronized pairs ( (</>;, <t>i+i) s.t. D(0i, <fii+i) ^ a) is tracked in a numerical 
experiment. In Fig 4, the x axis represents time, while y axis represents the count of neighbors that are 
not a-synchronized. From this figure we see that the trend of this system is achieving a-synchrony. In 
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Figure 4: Count of pairs in asynchrony with the change of time. The system parameters are a = 0.1, 
p = 0.3, e = 0.01, N = 200 



addition, a figure similar to [4] 's grey level denotation is made. In Fig 11 on page 15, x axis is the index 
of oscillators (that is to say, two pixels horizontally adjoining each other are neighboring oscillators i and 
i + 1), y axis represents time (time increases via the direction of up-down), and the greylevel of a pixel 
(i, t) on this figure depicts the phase fa of oscillator i at time t (the whiter, the nearer to 1, the blackest 
represents <f> = 0). This figure is a vivid record of the whole synchronization process. 

b. Testing time cost to achieve a-synchronization 

An experiment comparing the time used to achieve ^-synchronization with different system parameters 
is conducted. In the simulation, we tested systems with N oscillators, and repeated the test for 100 times 
on 100 randomly generated samples, for each parameter combination. We generated 100 initial condi- 
tion samples {Si,S 2 , • • • , Sioo}/ each S k = {(f>Q k \ fa[ h \ . . . , </>^_i}- For each combination of parameters 
such as a, p, e, an average time used to achive ^-synchronization of the 100 samples is calculated. See 
algorithm (1) for a more strict description of the algorithm. A result is showed in Fig 5 on the following 

Algorithm 1 Calculating the time cost to achieve a-synchronization with differnt p and e 
Randomly generate 100 initial samples {Si, S2, • ■ ■ , Sioo}. 

Generate the set of parameter values: p = {pi, P2, ■ ■ ■ , pi},e = {ei, £2, . . . , e m }. 

for all p e p do 

for all 6G? do 

timecost <— 0; 

for all 1 < is: 100 do 

timecost += (the time used to achieve a-synchronization by Si under p, e) 

end for 

_ . , timecost, 

Draw a point (p, e, — — Jon the p — e — t graph, 
end for 
end for 



page. 

c. Phenomena of asynchrony 

It is not always that, with any parameter, the system can eventually achieve a-synchronization. Similar 
to algorithm 1, a program to count the configurations that fail to achieve a-synchronization with a 
certain combination of parameters is designed. From Fig 6, it is easy to see that with a e [0.01, 0.1], 
asynchrony frequently happens. 

For a clearer observation, a 2-D plot with respect to a -count is drawn on Fig 7. Obviously, with 
a w 0.02 and a w 0.05, the chance of an asynchrony is large, while with a ^ 0.1 the experiment always 
end in synchrony. In purpose of trying to find a reason to these phenomena, an analysis of the whole 
section IV. is implemented. 





alpha 



Figure 6: A a — e— count graph that shows the 
count of configurations which don't achieve 
synchrony before t = 100. p = 0.3, N = 200, 
out of 100 samples. 
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Figure 7: Count of initial configurations that 
fail to achieve a-synchronization before t = 
100, out of 100 samples, p = 0.3, e = 0.01, 
N = 200. 



IV. Analysis of the Dynamics 



In purpose of explaining phenomena observed in Fig 7 on the page before, a step by step analysis is 
made in this section. First, we will review some useful notions introduced in previous works; then a 
new probability-based approach will be presented to find a cause of Fig 7 on the preceding page. 



a. Firing map 



Similar to [3], first we check the dynamics of a system of two oscillators: A, B. Now Let a row vector 
(4>i, (j)2) to represent the state of two oscillators. Without loss of generality, we can always "shift" the 
phases of A, B into a form of (a, 0), note that in this representation we change the order of components 
to make sure 4> > a. After a FD process, the state of this pair becomes (a,<f>') (a change of order from 
(4>', a)). We define function h : [0, 1) — ► [0, 1) to describe the mapping <j>' . Without considering the 
bounds, h may simply be 

h{4>) = {l + p)-{l + a-4>)+e (1) 
Now consider the restrictions on h(4>): 

1. when <f> G [0, a), there's no definition on h since we suppose > a. But for convenience we expand 
the domain of h to [0,1) 

2. when 4> € [a, 2a), since they achieve a-synchronization, there's no effect on A from B's firing. 

3. when <fi and a are "farther" away from a, but still too "close", B's firing will pull A's phase to reach 
1 and cause A's immediate firing. We can calculate this bound S a : 

h(5 a ) = 1 (1 + p) ■ (1 + a - 6 a ) + e = 1 <^> 

1 - e 



1 + a 



Based on the discussion above, we define the firing map 



1 + p 



(2) 



h(<f>) = { 



2a if < <j) < 2a 
if 2a ^ <f> ^ S a 

(1 + p) + </>) + £ if (5 a <</><! 



(3) 



b. Jumping length 

To inspect A's "length" of phase that was "pulled" forward by B's firing, we introduce the concept of 
Jumping Length. When A responds to B's firing at time t, A's jumping length is E>(</4 + , 4> l A ). Furthermore, 
on a standard two-oscillator configuration (a, cf>), the jumping length that A will have due to B's firing is 
defined as 

if 4> e [0, 2a) 

J(<f>) = I - 2a if e [2a, 5 a ] (4) 

p-{l + a-4>)+c if e (5 a ,l) 



And for a non-standard configuration (0a, 0b), where (0a < 0b)/ the jumping length of A that will be 
caused by B: J* (0a <— 0b) is defined as 

J* (0a <- 0b) = J (0b - 0a + a) (5) 
c. Stability of a-synchronization 

In a system of two oscillators, once the system achieves a-synchronization, the synchrony will never be 
lost as the oscillators don't respond to each other's pulse all the time regarding to the rule delay in the 
model. But in the system of more than 2 oscillators, each one is influenced by its two neighbors. Will a 
pair of oscillators that achieves a-synchronization lose synchrony due to others' coupling? Here we can 
show that if a pair is a-synchronized once, they will restore their synchrony any time one of them reacts 
to the other one's firing. 

Theorem 1. 0j and i+ i (0j) are neighboring oscillators that achieve a-synchronization at one time: 
3t : B(0j(to)> 0i+i(*o)) ^ a - Then for any t* > to at which ior i+ 1 fires, we will have 
D(0i(t* + a), i+ i(t* + a)) < a again. 

Proof. Without loss of generality, we suppose 4>i+\ = <f>i + d, ^ d ^ a, and k G {i, i + 1} the oscillator 
that will first fire (note that i may "catch up" with i + 1). Let t* be the time of fc's firing. 

Case 1 — Neither i or i + 1 receives any pulse in [to — a, t*\ : 

Obviously without the influence of pulse outside 

the pair i and i + these two oscillators will remain 
in synchrony: 

©(&(*)), &+i(t)) ^a,Vt G [t ,t*+a] 

(See the figure on the right: the height of bars repre- 
sents phase 0, the two bars in the middle are i and 
i + The leftmost and rightmost bars will not fire 
before i + l's firing.) 
Case 2 — i + 1 receives a pulse and i doesn't receive one before <*: 

Suppose i + 2 fires at t±, then i + 1 reacts to its pulse at £2 = ti + ol 

(i){l + P )-<P% l +e^l: 

As the phase of i + 1 reaches 1, i + 1 fires immediately and resets to (by the way we have 

t* = h) 

After a time length of a, i reacts to the firing at t 3 = t* + a. 

V 0i + i = 0i + d, d < a, 

■•■ 0* 3 = 0f + a = 0Ci - d + " > 0m 

=+ 4 = (1 + P) ■ 0? + e > (1 + p) • 0£i + e = (1 + P) ■ m +OI 
That is to say, 0j is reset to at t* + a, 
=^D(0i,0 i+ i)=D(O,a) = a, 

which indicates i and i + 1 are a-synchronized again. 




(ii)(l + p) ■ $ +1 + e < 1: 

After a time length I = 1 + a — 4> l 2 +1 , i will respond to i + l's firing. 
Let t 3 = t* + a ( also = t2 + Z), then 

(l + p).^3 +£ = (i + p ).^+l) + e 

= (l + p)-(<t>f +1 -d + a + l-cf>f +1 ) + e 

> + + (6) 

= i 

Hence </>j is reset to at t£, which implies E) (</>», 0»+i) = D(0, a) = a after i responds to the 
pulse from i + 1 

==> the ^-synchronization remains. 
So in case 2, the a-synchronization will restore. 

Case 3 — i receives a fire before t* 

This case is divided into two sub-cases with respect to k (whether i "catches up" with i + 1 or not). 

(i) k = i + 1 

i + 1 will be the first that fires and i will reacts to its firing at t 3 = t* + a 

Imagine there is a "shadow" oscillator i, which has a phase equal to i at to: 4>i ° = 4^ > but i 

doesn't respond to i — l's firing during [to,t* + a]. 

Then i 's case is discussed in case 2. 

~t+ 

=► 0i > 1. 

- — • ^ — 

And obviously we have < <^ 3 

& ><f>i > h 

which is a proof a i and i + l's restoration of a-synchronization. 

(ii) fc = i 

If z + 1 doesn't receive a fire from i + 2 before £3 = t* + a, 

Similar to Case 2(ii), suppose i — 1 fires at ti, then i reacts to its pulse at ti = t\ + a 
After a time length I = 1 + a — 4>? ,i + l will respond to i's firing. 
No matter i + 1 reacts to i + 2 's pulse before t* or not, we have 

(l + p)-^ 1 +e > (l + p)-(^ +1 +l) + e 

> (l+p)-(<f>f +d + a + l-4) + e 

> {l + p)-{4>f) + e+(l-4) 

= 4 + {i-4) 
= 1 

ThatisB(^,<^=B(a,0) = a. 
.'. the a-synchronization remains if i receives a fire before t* , no matter i + 1 receives or not. 



So we have proven that an a-synchronization will restore at the time when one of the once-synchronized 
oscillators reacts to its synchronized partner's pulse. □ 



d. Estimation on Dynamics 

The dynamics of a system on a bidirectional ring with delay is complex: each oscillator in it responds 
to both side of neighbors, while each neighbor is influenced by a farther neighbor. It's not easy to trace 
the evolution of the system based on an exact event-tracing approach like in [1], or [4] as the order of 
firing will not change due to unidirectional connection and non-delay in their cases. So, in this model, 
we need to use a new probability based approach to study the dynamics of the system. 

d.l Neighborhood of three 

Now consider a three-oscillators-composed neighborhood (B, C, D), in which B and C have not achieved 
a-synchronization (because the purpose of this analysis is on the synchronization process of B and C). 
Let a standard configuration (a, 0, 0), (0 ^ 2a)be used to describe the starting state of this neighborhood 
at to. For a given 0, if is in some certain area(which will be discussed in the next paragraph), C will 
have a coupling effect on B in current cycle. That is after a time of (1 + a — 0), is added up with 
J*{(j)^4>) — i.e. the state becomes (a + (1 + a - 0), + (1 + a - 0) + J*(0^0),a) at t + (1 + a - 0). 

d.2 Domain of 

For a given 0, there is an interval which must be in, in order to affect in C's current cycle. is 
defined piecewise: 

(0 + a, 1) U [0, a) if e (2a, 1 - a) 

(a + (j> - l,a) if (j) e (1 - a, 1) 

To unite the definition, we extend <f> 's domain from [0,1) and treat 1 + the same as 0. In this way 
becomes 

D = (0 + a, 1 + a) 

d.3 Reducing three to two 

With respect to D's influence on C, there is no difference whether C's "adding up" J r *(0^0) happens 
after a while or immediately. So, in purpose of analyzing the dynamics of (B, C) in one cycle, we can 
treat D's coupling as if it happens immediately at to. That is, for this neighborhood (a, 0, 0), we can 
depict its state with one less opponent: (a, + J* (0«_0)) , by making D's pulse into effect immediately 
and discard it from this neighborhood for its effect in this cycle is already done. So the representation of 
this neighborhood (a, 0, 0) becomes (a, + J"*(0^0)). 

d.4 Random factor 



By putting the neighborhood (B, C, D) back to the ring of oscillators, we may notice that B is also influ- 
enced by its left neighbor A, while D is affected by its right neighbor E. If we go on inspecting A and 




Figure 8: A demonstration of replacing 0d with the mathematical expectation of influence that it may 
cause on C. 

E, we have to face a very long chain of cause-and-effects. So, now we use a probability based approach 
to analyze the possible behavior of the system. That is, when analyzing the behavior of B and C, we 
consider A and D as random factors — their phases uniformly distribute on [0,1). 



d.5 Expectation of influence 

For a standard configuration of B, C: (a, </>), C is possible to be influenced by D's firing before its reaching 
0. So we consider D's phase <j> a random variable uniformly distributed on [0, 1). When <f> 6 T>^, C will 
react to D's firing before its own resetting. This response causes C's phase to have a jumping forward of 
length J*{<j)^(f). So the mathematical expectation of C's jumping length caused by random variable cf> 
is 



1 



J* 



1 - 

J*(<t>^<t>)d<f> (7) 



So in the neighborhood consisted by B, C, D (a, (f>, (f>), we regard D as a random factor and replace it by 
its "average" affection E(J*{(f>^ip)) on C. Thus we can analyze B,C's dynamics in a cycle by inspecting 
(a, 4> + E{J*{(f)^(f)))). This replacement is demonstrated in Fig 8. 



To calculate + E{J*(<j)+_<j))) using (7) and (5), we have 

J*U^$) 



l+a 

J(4>-<f> + a)i 
J(s) ds 



(s ~ 2a) ds + / (p ■ (1 + a - s) + e) ds if 1 + 2a ~ <\> > S a 



2a 
l+2a- 



(s - 2a) ds 



That is (see Fig 9 for a plot) 

if <h < 1 + 2a - S a : 



if 1 + 2a - 6 s: So 



■ + E{J*{4>^4>)) = < f>+l.(8 a - 2a) 2 + ~ (1 + 2a - 5 a - cf>)[2e + (l-6 a + cj>)p] 



if <j> > 1 + 2a - 5 a : 



c^ + E{J*{^cj>)) 



figure of E{J*(<p<-<!> )), alpha=0.1, rho=0.3, epsilon=0.01 



Figure 9: A plot of 
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o* j E(J*(<fn-<l>)) 



d.6 Absorption of a-synchronized pairs 

Suppose there is a pair of a-synchronized oscillators A and B, whose states can be depicted in a standard 
form (<fio, a) where 4>q = a—d, ^ d ^ a (we can always make <po > a with symmetry of the bidirectional 
ring). Oscillator C (</>i) is B's right side neighbor. If ID(0c, 4>b) > a, we will find out that B and C are 
likely to achieve a-synchrony, which looks like C is absorbed into an a-synchronized pair near it. 

Because A and B have already achieved a-synchrony according to Theorem 1, this pair will restore 
their a-synchrony everytime they reacts to each other's pulse. In this estimation, we ignore the small 
chance that A may catch up with B, so B only reacts to C's pulse. Follow the mathematical expectation 
based approach, we may estimate B and C's state as (a, <f)+E(J* ((/><_ cf>))). After B's reacting to C's firing, 



the state becomes (h(4> + E(J* (</>«_0))), a). Adding an expectation of D's influence on C again, the state 
goes to 

(htf + E(J*(<t>U))), a + E(J*(aU))) 

which is equivalent to 

(htf + E(J*(^))) - E(J*(aU)), a) 

in the next round it becomes 

a,h(h (</> + E(J*(<f>^j>))) - E(J*(a^$)) 

Note that J*((t> <— 4>) defined in (7) has a prerequisite of 4> > 2a, so we have to recalculate: 

E(J*(a^4>)) 

= E{j{4>)) 
= I J{s)ds 

= / (s-2a)ds + (p ■ (1 + a - s) + e) ds 

Jla JS a 

= X - (S a - 2a) 2 + (1 - S a ) [p (1 + a) + e] - \ (1 - S a ) 2 

We define the probability firing map: 

M{<j>)=h(h{j> + E{J*{^4>)))-E{J*{a^)) S j (8) 

Then similar to [1] and [3], the synchronization of B and C can be studied by analyzing the iteration of 
H. That is, whether C will be absorbed to an a-synchronized pair (A,B) depends on whether leN s.t. 

H( fe )(0) = 0, here is an iteration of H: = H(H( fe - 1 )) 

d.7 Explanation of asynchrony phenomena with iteration of H 

In purpose of analyzing the properties of H, its graphs are plotted with respect to different param- 
eters. We noticed the change of the structure of H 's graph with different a. Also, in the experi- 
ments mentioned in previous section, we observed that with a ^ 0.1, systems often fail to achieve 
^-synchronization. With the graphs of H with different a, such phenomenons are easier to explain. 
From Fig. 10 on page 17, we can see that with a < 0.1, the H 's graph is primarily consistied of a skew 
part and two flat parts, defined as = {x|H(.x) = 0}. Obviously, if x G N^, B and C will achieve 
synchrony in the next round. Even if x ^ TVh, as the slope of the skew part [0, l]\ATn is greater than 1, 
the value of M k (x) will eventually "fall" in to A^ H , which results in M k+1 (x) = 0. 

As the value of a increases, the graph of H(x) gradually becomes not so good looking. For a = 0.02 in 
Fig 10 on page 17, we can find "pits" on the graph, which prohibits the iteration to reach M k (x) = 0. For 
example, the flat part [0.1, 0.15] will cause every initial values that in it to reach a fixed point x* w 0.11, 
i.e. Vx e [0,0.15], H(x) = x*, where H(x*) — x* ^ 0. This is an explanation of the happening of 
phenonemons of asynchrony in a = 0.02 observed in previous experiments in fig 7 on page 6. 



)) 



Also on the graph of a = 0.05, the same patterns are spotted: H([0, 0.15]) « 0.1 e [0,0.15], which 
creates a fixed point of H that attracts initial values in [0, 0.15]. 

However, with a = 0.2, the graph of H in Fig (10) is in another pattern. A rough observation on this 
graph tells: for turning points x\ ~ 0.25, x 2 « 0.4, x 3 w 0.75, x 4 w 0.9 of y — H(x), H([0,xi]) > xi, 
H([xi, x 2 ]) > x 2 , H'([x2, X3]) > 1 and H([x 3 , x 4 ]) = 0, H([x 4 , 1]) e [0, xi], which means 



[x 4 ,l] ^ [0, Xl ] 



JH 



[0,Xl] ee^- [xi,x 2 ] 



[xi,x 2 ] ^ [x 2 ,x 3 ] 



[x 2 ,x 3 ] 



So, for an arbitrary initial point (j> e [0, 1], with iteration of <j) k+ i = H(0 fc ), we will inevitably have 
4> n = for some n e N. 



V. Conclusions 

In this paper we have proven that in each pair that once achieved 
a-synchronization will restore their synchrony no matter how 
their neighbors acts. Then with comparison between the numer- 
ical simulation results and the properties of H, we have success- 
fully established a cause-and-effect link between the synchro- 
nization behavior of the systems and their probability firing map 
H. With this application of H, it is reasonable to believe that the 
probability based approach presented in section (IV.) reaches the 
point of the problem. 

VI. Further works 

On Fig 5 on page 6, a non-linear change of timecost with respect 
to a is clearly seen. But this phenomenon is not yet studied in 
this paper. Also, in section (IV.), we ignored a lot of affecting 
factors that may attribute to the result of the system's evolution. 
And what is more, the reason of forming of patterns in H should 
be searched with analytic studies on E(J*(4> <— <fi)) Besides, we 
would like to inspect synchronization phenomena when a central 
"host" oscillator is added, which has connection to all, or part of 
the oscillators on the circumference. We hope to work on those 
topics in further studies. 
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Figure 11: This figure 

demonstrates an a- 
synchronization process 

of 100 oscillators, x axis 
represents the index of 

oscillators, y axis rep- 
resents time, grey level 
shows the phase of the 
oscillator (blackest: <p = 0, 
whitest: 0=1). Inspired 
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Figure 10: Graph of H(</>) with differnet a, when p — 0.3, e = 0.01 



